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Delocalized Operator Expansion 
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A generalization of Wilson's local OPE for the short-distance expansion of Euclidean current correlators, called 
delocalized operator expansion (DOE), which has been proposed recently, is discussed. The DOE has better 
convergence properties than the OPE and can account for non-local non-perturbative QCD effects. 



1. INTRODUCTION - WILSON'S OPE 

The Wilson OPE is one of the standard tools in 
modern hadronic physics. The OPE provides the 
framework for systematically separating short- 
distance contributions (x <C A -1 )^ from long- 
distance contributions (x ~ A -1 ). For illustra- 
tion consider the large-momentum expansion of 
the correlator of a gauge-invariant QCD current 
Six) [1], 

Q 2 EX>™W 2 )<°w(o)) (1) 

N=0 1=1 

where N runs over the dimension of the local 
and gauge-invariant composite operators Ojvi(0), 
and I labels operators of the same dimensions. 
The long-distance fluctuations are encoded in the 
matrix elements (Ojv/(0)), also called conden- 
sates, and the short-distance contributions are 
contained in the Wilson coefficients cni{Q 2 )- One 
has the scaling (Ojvz(O)) ~ A N and cni{Q 2 ) ~ 
Q~ N . The N = term is entirely perturbative 
and higher order terms in the series are of or- 
der (A/Q) . For the case Q 3> A the expan- 
sion should be reasonably well-behaved, but it is 
known to be asymptotic. There are situations 
where the OPE cannot be applied or where it can- 
not give definite answers. When Q approaches A, 
the series breaks down. The OPE can also not an- 
swer questions related to the asymptotics of the 
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series due to the truncation of the series. In this 
talk I discuss a generalization of the local OPE 
that has been proposed in Ref. [2] and has been 
called "Delocalized Operator Expansion" (DOE). 
The DOE is an attempt to combine the advan- 
tages of the OPE with methods that might be 
useful to resolve issues that cannot be tackled eas- 
ily in the framework of the OPE. 

2. BASIC FRAMEWORK 

Given the exact knowledge of (nonlocal) gauge- 
invariant QCD vacuum correlators and assum- 
ing for now that perturbative propagation can be 
factorized unambiguously from nonperturbative 
fluctuations, one may view the OPE as the mul- 
tipole expansion of the perturbative part. To see 
this, let us consider the simplified case of a non- 
perturbative and nonlocal dimension-4 structure 
g(x) corresponding to a slowly varying 2-point 
correlator falling off at distance A g ~ A -1 in its 
Euclidean space-time argument x. In addition, 
we consider a perturbative short-distance func- 
tion f(x), that probes the vacuum at distance 
x ~ A f ~ Q~ l . For Q > A, f(x) is strongly 
peaked compared to g(x) (see Fig. 1). For illus- 
tration I consider the problem in one dimension. 
Then the chain of local power corrections to the 
perturbative result is obtained from the expres- 
sion 

jr oo dxf(x)g(x)=:(f,g), (2) 

which is the bilinear form in the dual space 
spanned by strongly peaked functions such as / 
and slowly varying functions such as g. The gen- 
eralization to more than one dimension, which 
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Figure 1. Schematic picture of the short-distance 
function f(x) and the 2-point-correlator g(x) illus- 
trating the scale hierarchy A//A 9 ~ A/Q <C 1. 

also allows the treatment of non-perturbative n- 
point correlators, is straightforward. In the OPE 
f(x) is expanded in analogy to the multipole ex- 
pansion of a localized charge distribution in elec- 
trostatics. This can be formulated by introducing 
the dual space basis functions (n = 0, . . . , oo) 



!„(i) = ^iW(i), e n (x) 



(3) 



with the orthonormality relation (e n , e m ) = 5 nm . 
The multipole expansion of f(x) then reads 



f( x ) = T,n=o(f^n)e n (x) 

which leads to 



(4) 



„ OO OO 

dxf(x)g(x) = ^2{f,e n )(e n ,g) =: ^f n g n . (5) 

•* n=0 n=0 

The /„'s are the Wilson coefficients and the g n 
the matrix elements of operators obtained from 
locally expanding the field content in the correla- 
tor g(x) for small x. In momentum space repre- 
sentation the Wilson coefficients and matrix ele- 
ments have the generic form 

fn = /»k=0 , 9n = J f ^ §(k) , (6) 

where / and g are the Fourier transforms of / and 
g, respectively. 

Briefly switching to 4 dimensions, a simple ex- 
ample for the function g(x) is the gauge invariant 
field strength correlator [3] 

9llVK x{x) = Tr(. 9 2 G^(x)5(x, 0)G kA (0)S(0, a;)), (7) 

where 

S{x,y)=V exp{ Jff /; dz^(z)}. (8) 



It generates the following chain of vacuum expec- 
tation values (VEV's) involving local operators of 
increasing dimension 



Tr(g 2 G^G KX ), Tr(g 2 G^D p G KX ), 
Tr(g 2 G^D T D p G KX ) 



(9) 



In the limit, where R = Ag/Af — > oo, the 
VEV with the lowest dimension, usually called 
the gluon condensate, dominates, and the con- 
tributions of higher dimensional VEV's are sup- 
pressed by higher powers of A/Q. VEV's with 
odd numbers of covariant derivatives do not con- 
tribute due to parity and time-reversal invariance. 
Pcrturbativcly, one can define the gluon conden- 
sate in such a way that its anomalous dimension 
vanishes to all orders in a s . 

Improved convergence properties may be 
achieved using a multipole expansion of f(x) 
based on functions of width Af w Q^ 1 instead of 
the infinitely narrow (5-functions. This is the basic 
idea in the construction of the DOE. At this point 
I would like to note that there exist a number 
of phenomenological studies where the non-local 
expression in Eq. (2) has been analyzed directly 
without any expansion for cases where the OPE 
did not have a good convergence behavior, see 
e.g. Refs. [4]. This approach has the feature that 
it requires a model-dependent ansatz for the cor- 
relation function g(x) and that the computation 
of f(x) can become cumbersome, particularly at 
higher loop level. The DOE has been constructed 
with the aim to provide an alternative formal- 
ism to describe non-local effects. As we will see 
later the DOE simplifies numerical predictions in 
a given model for g(x). However, I will also show 
that the DOE allows to extract non-local non- 
perturbative information on the QCD vacuum in 
a model-independent way. 

Let me continue with the construction of the 
DOE. For Cartesian coordinates the dual space 
basis 



W) = ^H n (Qx)e 

e nW — (2[))» ! 



(10) 



with (e^,e^) = 5 mn , the H n being the Hermite 
polynomials, is well suited. Of course this choice 
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is not unique, but it fixes a scheme, which can be 
unambiguously related to other possible schemes 
that can be used. The have a width of order 
and for U — > co one finds — > e„, — ► e„. 
The parameter f2 is called resolution scale. In this 
basis Eq. (2) can be written as 

£~=0(/>3?)(«£,S) = : E^o/™(")5n(^), (11) 

where the fi-dependent short-distance coefficients 
and matrix elements have the form 

9n m=j£&£e-^g(k), (12) 

in momentum space representation. The series in 
Eq. (11) has better convergence properties, if is 
of order Q rather than being equal to oo, because 
the first term in the delocalized multipole expan- 
sion generally provides a better approximation of 
the actual form of f(x) than the local expansion. 
The relation between basis functions for resolu- 
tion parameters and tt' reads 

^W=E™=0«n m (O,O')^'W 

#(aO=£~=o<W fi '> fi ) (13) 
where 

/ n—m \ 

n\ /SI 2 — Sl' 2 \ v 2 ^ 

a nm (n, n') = ^(a^m), ^4^2 n , 2 J ( 14 ) 

if n — m > and even, and a nm = otherwise. 
The transformations a nm (Cl,Q') form a group. 
One can relate the short-distance coefficients and 
matrix elements for different resolution parame- 
ters to each other. Here, I only want to discuss 
these relations for finite £1 and CI' = oo: 

/o(fi) = fo , (15) 

/2(0) = /2 - 5?f/o , • • • j 

f fr)\ _ V^[™/ 2 ] rt! I jr 

Jn\ il -) — L-ii—o (n-2i)\i\\ 4fi 2 / J n-2n 

5n(«)=£^^(^) l 5n + 2,. (16) 

The coefficient / n (f2) can be expressed in terms 
of a /mite linear combination of the local Wil- 
son coefficients /,(oo) for i < n. The short- 
distance coefficient of the leading power correc- 
tion is f2-indcpcndcnt. The ^-dependent matrix 



elements <7 n (f2) are related to an infinite sum of 
local matrix elements with additional covariant 
derivatives. These properties exist for all dual 
space bases constructed from orthogonal polyno- 
mials. Note that relations (15) can also be cm- 
ployed if the separation between long- and short- 
distance contributions is factorization scale de- 
pendent, and one can therefore use these relations 
as the formal definition of the terms in the DOE. 

The DOE of quantities such as in Eq. (1) has 
the same parametric counting in powers ofA/Q as 
the OPE, if is not chosen parametrically smaller 
than Q. Consider that g n ~ A™ and /„ ~ Q~ n in 
the OPE, then we have 

s„(fi)~A» £.(£)' ~A», 
/^^-"E^f-Q-", (17) 

as long as > Q. However, one expects that the 
actual size of the term f n (£i)g n (£Y) in the DOE for 
fl ~ Q has an additional numerical suppression 
by powers of a small number. 

3. HEAVY QUARKONIUM GROUND 
STATE ENERGY 

In heavy quarkonium systems the relevant 
physical scales, mass m, momentum p, energy E 
and A have the hierarchy 

m > p ~ mv > E ~ mv 2 > A, (18) 

where v <C 1 is the quark velocity. Thus the spa- 
tial size ~ (wifj) -1 is much smaller than the typi- 
cal dynamical time scale ~ (mv 2 ) -1 . In this sec- 
tion I demonstrate the DOE in a toy-model com- 
putation for the nonperturbative corrections to 
the l 3 5i ground state for the expansion in A/E. 
My intention is not to carry out a phenomeno- 
logical study, but to show how well the series in 
the DOE behaves in comparison to the series in 
OPE and how well the DOE approximates the 
exact model result. The ratios of m, p and E 
are treated at leading order in the local expan- 
sion. This means that the perturbative dynam- 
ics is described by the nonrelativistic two-body 
Schrodinger equation and that the interaction 
with the nonperturbative vacuum is accounted for 
by two insertions of the local xE dipole operator, 
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Figure 2. The perturbative short-distance function 
f(t) for m — 5 GeV and a s = 0.39 (solid line) and the 
leading term in the delocalized multipole expansion of 
f(t), [J dt'f(t')e%(t')] for Q = k 2 /m (dashed line). 

E being the chromoelcctric field. [5] The chain of 
VEV's of the two gluon operator with increas- 
ing numbers of covariant derivatives times pow- 
ers of quark-antiquark octet propagators [5], is 
treated in the DOE. In this model the interaction 
with the vacuum fluctuations only depends on the 
temporal distance (in Euclidean space) of two in- 
sertions of xE dipole operators and the nonper- 
turbativc correction to the ground state energy 
reads (k = |ma s ) 

E np = S- X dtf{t)g{t), 
.f(t)=/^ e -^/d 3 x/d 3 y^(x)(xy) 

xG o (x,y,-£-q H(y), (19) 

where f(t) is the short-distance function that de- 
pends on ground state Coulomb wave function 
</> and the octet Green-function G [5] and g(t) 
is the gluon field strength correlator for x = 
(t, 0,0,0). In Fig. 2 the function f(t) (solid line) 
is displayed for m — 5 GeV and a s = 0.39. The 
characteristic width of / is of order the energy 
k 2 /m ~ (ma 2 )^ 1 ~ (m» 2 ) -1 . For the nonper- 
turbative gluonic field strength correlator I use a 
lattice-inspired [6] model 

g (t) = 12 Aq exp(- ^t 2 + X\/X A + 1) 

A = 0.04 GeV 4 , A A X = 0.7 GeV , (20) 

with a large-time behavior ~ e ^ t / ly A j n ^j^ig 
model the value of the gluon condensate is 
(— G" } = ^4* = 0.024 GeV 4 . In Tab. 1 
the exact result and the first four terms of the 
resolution-dependent expansion of E np are shown 
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m 




(GeV) 


E np 
(MeV) 


n 


fnQn 

(MeV) 


fnQn 

(MeV) 


0.39 


0.338 


24.8 





38.6 


24.2 








2 


-65.7 


-3.9 








4 


832.7 


12.1 








6 


-35048.0 


-43.1 


0.15 


1.750 


0.245 





0.258 


0.249 








2 


-0.016 


-0.005 








4 


0.008 


0.003 








6 


-0.012 


-0.003 



Table 1 

Nonperturbative corrections to the heavy quarko- 
nium ground state level as described in the text. 



for the quark masses m = 5 (upper part) and 
175 GeV (lower part) and for ft = oo and ft = 
k 2 /m. The strong coupling has been fixed by 
the relation a s = a s (k). The numerical values 
of f n g n have been determined from Eqs. (6,12). 
The series are all asymptotic. The local expan- 
sion is badly behaved for m — 5 GeV because 
k 2 /m < , and basically meaningless. For 
m = 175 GeV, where k 2 /m > , the local ex- 
pansion is good. For Q = k 2 /m, however, the se- 
ries is much better behaved for all quark masses. 
The size of the order n term is about a factor 
2~™ smaller than the order n term in the local 
expansion. One also observes that even in the 
case k 2 /m < A^ 1 the leading term in the delo- 
calized expansion for ft = k 2 /m agrees with the 
exact result within a few percent. This feature is 
a general property of the delocalized expansion, 
and should apply to any quantity for which the 
local expansion in the ratio of two scales breaks 
down because the ratio is not sufficiently small. 

4. RUNNING GLUON CONDENSATE 
FROM C HARMONIUM SUM RULES 

The Sl-depcndcnt matrix elements are either 
determined from experimental data or from lat- 
tice measurements. In the following the in- 
dependent ("running") gluon condensate is ex- 
tracted from charmonium sum rules which are 
based on moments [7] 

M n = M-m*) ni i c (Q 2 )\Qi=o (21) 
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Figure 3. The running gluon condensate as a func- 
tion of n when extracted from the ratio of charmo- 
nium moments r n . 



of the correlator of two charm quark vector cur- 
rents j M = C7 M c. The moments can be determined 
theoretically in an expansion of the form [8] 

M n = M° n {\ + [p. corr.] + 6^ (g 2 G 2 ) + . . .} ,(22) 

while the experimental moments are obtained 
from a dispersion integral over the cc cross sec- 
tion in e + e~ annihilation. The Wilson coefficient 
of the gluon condensate, S^n \ is ^-independent. 
We consider the ratio [7] r n = ^ n - and extract 
the gluon condensate as a function of n. Since 
the relevant short-distance scale for the moment 
Ai n is of order m c /n, a proper choice for the res- 
olution scale is O = 2m c /n. Thus the depen- 
dence of the gluon condensate on n can be inter- 
preted as the dependence on ft. The results for 
the gluon condensate as a function of n is shown 
in Fig. 3 for n < 8 and m c (rn c ) = 1.23 (white tri- 
angles), 1.24 (black stars), 1.25 (white squares) 
and 1.26 GeV (black triangles). The area between 
the upper and lower symbols indicates the exper- 
imental and theoretical uncertainties. The run- 
ning gluon condensate appears to be a decreasing 
function of n. The solid thick line is (— G 2 )(f2) 
obtained from the form of the gluon held strength 
correlator suggested from lattice computations [6] 
for Q = (2.5 GeV)/n and the vacuum correlation 
length A^ 1 = 0.7 GeV. The qualitative agreement 
is encouraging but, the uncertainties of our ex- 
traction are still quite large. (See Ref. [2] for an 
analysis based on hadronic r decay data.) 



5. CONCLUSIONS 

In this talk I discussed a generalization of the 
Wilson OPE based on nonlocal projections of 
gauge invariant correlation functions in a delo- 
calizcd version of the multipolc expansion for the 
perturbatively calculable coefficient functions. [2] 
This "Delocalized Operator Expansion" (DOE) 
depends on an additional parameter ft, called 
"resolution scale" which adjusts the width of the 
projection functions used for the expansion. The 
DOE has the same power counting as the OPE, 
but in general better convergence properties than 
the OPE. The DOE short-distance coefficients 
can be determined from the OPE Wilson coeffi- 
cients at the same order in the expansion, whereas 
the DOE matrix elements correspond to an infi- 
nite sum of OPE matrix elements with additional 
covariant derivatives. I believe that the DOE can 
serve as a useful tool for situations where the OPE 
cannot be applied. 
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